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Fabrizio Altarelli, 1,2, EI Alfredo Braunstein, 1,3,2 Luca Dali’Asta, 1,2 Caterina De Bacco, 4, @ and Silvio Franz 4 

l DISAT and Center for Computational Sciences, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy 
2 Collegio Carlo Alberto, Via Real Collegio 30, 10024 Moncalieri, Italy 
3 Human Genetics Foundation, Via Nizza 52, 10126 Torino, Italy 
*LPTMS, Centre National de la Recherche Scientifique et Universite Paris-Sud 11, 91405 Orsay Cedex, France. 

We present a message-passing algorithm to solve the edge disjoint path problem (EDP) on graphs incorporat¬ 
ing under a unique framework both traffic optimization and path length minimization. The min-sum equations 
for this problem present an exponential computational cost in the number of paths. To overcome this obstacle we 
propose an efficient implementation by mapping the equations onto a weighted combinatorial matching prob¬ 
lem over an auxiliary graph. We perform extensive numerical simulations on random graphs of various types 
to test the performance both in terms of path length minimization and maximization of the number of accom¬ 
modated paths. In addition, we test the performance on benchmark instances on various graphs by comparison 
with state-of-the-art algorithms and results found in the literature. Our message-passing algorithm always out¬ 
performs the others in terms of the number of accommodated paths when considering non trivial instances 
(otherwise it gives the same trivial results). Remarkably, the largest improvement in performance with respect 
to the other methods employed is found in the case of benchmarks with meshes, where the validity hypothesis 
behind message-passing is expected to worsen. In these cases, even though the exact message-passing equations 
do not converge, by introducing a reinforcement parameter to force convergence towards a sub optimal solution, 
we were able to always outperform the other algorithms with a peak of 27% performance improvement in terms 
of accommodated paths. On random graphs, we numerically observe two separated regimes: one in which all 
paths can be accomodated and one in which this is not possible. We also investigate the behaviour of both the 
number of paths to be accomodated and their minimum total length. 
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FIG. 1: An instance of the EDP problem over a 3-regular random graph of V = 20 and M - 6: examples of solutions of the 
unconstrained (left) and optimal (right) EDP problem are displayed. In the latter, the purple communication is redirected along 
a longer path to avoid edge-overlap. The yellow one has two shortest paths of equal length (degeneracy) in the unconstrained 
case, but once EDP is enforced this degeneracy is broken and only one of the two is optimal (right). 


I. INTRODUCTION 


The optimization of routing and connection requests is one of the main problems faced in traffic engineering and communi¬ 
cation networks Cl- The need to deliver Quality of service (QoS) (2 b 3| performances, when transmitting data over a network 
subject to overload and failures, requires both efficient traffic management and resource optimization. 

Some aspects of these problems can be formalized using the edge-disjoint path (EDP) problem. This is a constrained opti¬ 
mization problem that is defined as follows. For a given network and a set of communication requests among pairs of users, 
the EDP consists in finding the maximum number of communications that can be accommodated at the same time, under the 
constraint that different paths cannot overlap on edges. Moreover, the additional requirement of minimization of the total path 
length can be considered. Apart from a purely theoretical interest @1, the EDP finds a wide range of applications: in very-large- 
scale-integration (VLSI) design, in admission control and virtual circuit routing and in all-optical networks. In VLSI design it is 
required to route wires on a circuit avoiding overlaps, along with minimizing the length of the wires nil. In admission control 
and virtual circuit routing 1749] one needs to reserve in advance a given path for each communication request so that once the 
communication is established no interruption will occur. This has applications in real-time database servers, large-scale video 
servers ee, streaming data and bandwidth reservation in communication networks iTsUThll and in parallel supercomputers. 
All these applications require high quality data transmission and full bandwidth exploitation. Routing via edge disjoint paths 
allows for an efficient bandwidth allocation among users because overlap avoidance means full bandwidth exploitation by each 
single user. An area that has attracted particular attention in the last decade is communication transmission in all-optical net¬ 
works. Along an optical fiber different communications cannot be assigned the same wavelength to transmit data. Moreover, 
a unique wavelength must be assigned on all the edges contributing to the path assigned for a given communication. Routing 
communications under the above two requirements define the problem of routing and wavelength assignments (RWA) in this 
type of networks El- These two constraints suggest that a strategy that iteratively builds edge disjoint paths solutions could 
allow for a more efficient bandwidth management, namely by using an overall smaller number of wavelengths. This leaves 
available the remaining ones (according to the edge capacity) to be used either by new users entering the network or by allowing 
current users to exploit higher bandwidth. This strategy has indeed been applied using greedy [18] and genetic algorithms ill911 
with performances comparable to other methods based on integer linear programming, graph coloring or bin packing. 

The EDP is classified among Karp’s NP-hard combinatorial problems 1201 u ill . Defining the approximation ratio of a given 
algorithm as the ratio between the result obtained in term of cost/profit by the algorithm and the optimal one (or viceversa de¬ 
pending on what order gives the maximum ratio), the EDP problem is hard to approximate in the worst case; it has been proved 
that even an approximation with ratio 0(m^~ e ) is NP-Hard. The best known approximation ratio for the number of accomodated 
paths is (9(minj/A 3 , ^Jm}) |22[|23|] where n and m are the number of nodes and edges in the graph, respectively. Negative results 
on worst-case inapproximability did not stop progress on heuristc approaches. The problem has been studied intensely with a 
variety of classical techniques: heuristic greedy algorithms El 313 EH, elaborated strategies using bin packingEjl inte¬ 
ger/linear programming relaxations ll26l - l29tl . post-optimization pO l, Montecarlo local search Oltl . genetic algorithms im, 
particle swarm optimization |35il and ant colony optimization 13611 . among them. 

In this paper we propose a distributed algorithm to solve the EDP problem based on message-passing (MP) techniques (or 
cavity method) [!37{]. This method has been extensively employed to address problems in spin glass theory 113 81 - 14 (ill , combinatorial 
optimization 114111 and more recently in routing problems on networks ll42l - l45ll . The evaluation of the equations at the core of the 
MP technique requires, for each vertex i in the underlying graph, to solve a local combinatorial optimization problem, performing 
a minimum over a set which is exponentially large in the number of neighbors of i. We propose an efficient method to perform 
this calculation, by mapping it into a minimum-weight matching problem on a complete auxiliary graph with vertices in the set 
di of neighbors of i, that can be solved by classical algorithms I4al . With this construction, each iteration of the MP equations 
can be computed in a time which is polynomial in the number of graph edges (and linear in average for sparse random graphs). 
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The MP algorithm is tested on computer-generated instances of different classes of random graphs to study the scaling properties 
with the system size and to compare the performances against a greedy algorithm. We also considered the EDP problem on some 
benchmark instances found in the literature, for which we could compare the message-passing results with those obtained using 
other types of algorithms: greedy, ant colony optimization [|36j and Montecarlo local search Dill . 

The paper is organized as follows. In Section[II]we define the EDP optimization problem, for which we present the message¬ 
passing equations in Section Hill together with the mapping on a matching problem that simplifies their actual implementation. 
Section |Iy] reports the results of simulations on random graphs and the scaling of the relevant quantities with the system size, 
while the comparison between the performances of the message-passing algorithm and other methods is discussed in Section [V] 
Conclusions are given in SectionlvD 


II. THE EDGE-DISJOINT PATHS PROBLEM 


Given a network and a set of M communications requests between pairs of senders and receivers, the standard EDP consists 
in finding the maximum number of accommodated paths which are mutually edge disjoint. In the applications described in the 
Introduction, the length of the communication paths is a quantity that, directly or indirectly, affects the overall transmission 
performances, in terms of transmission delays, infrastructure cost and network robustness. We take into account this aspect by 
considering the Minimum Weight Edge Disjoint Paths (MWEDP) problem, a generalization of the EDP problem that combines 
in a unique framework both path length optimization and edge disjointness. An instance of the MWEDP problem is defined by a 
graph GOV, £), where *V denotes the set of nodes and £ is the set of edges, by an assignment of edge weights w, that we assume 
to be non-negative real numbers and by a set of M communication requests {(.S’ A ', between ordered pairs of nodes. 

We denote by n a path, i.e. a set of consecutive edges e e £, that connects a sender S 11 with the corresponding receiver R 1 '. The 
optimization problem consists in finding M pairwise edge-disjoint paths 7r |( while minimizing the total edge weight Z t t w(n tl ), 
where w(n M ) = Zeen„ Me)- 

The classical EDP problem could be trivially recovered by assigning zero weight to all edges in G(W, £) and a positive cost 
to each communication that is not accommodated. Alternatively, any solution of the MWEDP problem can be reinterpreted as a 
solution of the classical EDP problem by slightly modifying the original instance of the graph G('V, £) by introducing an extra 
edge between each pair (5^, 7?^) with sufficiently large cost, such that the algorithm could still always find a solution possibly 
using these expensive extra edges. By construction, the cost of each of these M extra edges should be larger than the maximum 
possible weight a single path can take. Then the solution of the classical EDP problem is obtained from any solution of the 
MWEDP problem by discarding the paths passing through the extra edges. In the present paper, we keep information about path 
length minimization by assigning unit weights (i.e. vv /; = 1, i(ij) e £) to the original edges of the graph G(*V, 3) and a fixed 
cost |£| + 1 to the extra edge added between each pair (S M , R>*). 

We introduce M-dimensional variables /,■>■ = (/?., ... ,lff) with entries /f. e {±1,0} representing the communication passing 
along an edge: 


I 


‘j 


1, if communication p passes from i to j, 
-1, if communication p passes from j to i, 

0, otherwise. 

7 


(1) 


We call these vectors currents as they must satisfy current conservation at each node i (Kirchhoff law): 

Y J 1 - j + A?=0, Vp = 1,... ,M, (2) 

jedi 

where we defined for each node i and each communication p a variable A 1 ' such that 


A? = 


if i - S M , 
if i = RZ 
otherwise. 


(3) 


The constraint of edge-disjointness specifies that for each edge ( ij ), at most one of I 1 '- is non-zero, therefore each vector Ijj can 
be parametrized by a variable taking 2 M + 1 different values. Notice that the set of variables {/qlqpe/; completely specifies the 
state of the network. In this multi-flow formalism, the MWEDP problem is a combinatorial optimization problem in which the 
global cost function C({/p}) = Zuj)eE w ijf(\\hj\\) depends additively on the total net current ||/ (; j| = |/f ; l along the edges, and 

the edge-disjointness is ensured by defining 


( 0 , 


miijW) = • 


i, 

+oo. 


if Poll = 0, 
if II4/H - 1, 
if ||/, 7 || > 1, 


(4) 
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where \I^\ = 0, 1 denotes the absolute value of I 1 ' . Thus configurations with more than one communication passing along an 
edge have infinite cost and, in the case of unit weights, the total cost C({/ |; }), if finite, represents exactly the total path length, i.e. 
the number of edges traversed by paths. 


III. THE MESSAGE-PASSING ALGORITHM 


On a tree, the optimization problem defined in SecQUcan be solved exactly by iteration using the following message-passing 
algorithm. Let us assume that G is a tree and consider the subtree G[ !y ] defined by the connected component of i in G \ (ij) (see 
Figure 0. We define Ejj(Ijj) to be the minimum cost C({/,,}) among current configurations that satisfy Kirchoff’s laws on all 
vertices of G[,/] given that we fix an input (or output) extra current l, f entering (or exiting) node i. Because of the absence of 
cycles, it is possible to write a recursive equation for £j y as a sum of cost contributions coming from neighbors of i in the subtree, 
plus the single cost contribution due to the current 7 ly passing along edge (ij). We call these quantities Ejj(Ijj) messages and they 
verify the min-sum recursion relation |[39i|: 


Eudij) 


min 

{Ikj)\constraint 



+ /(l|7 iy ll) 


(5) 


where constraint is the Kirchhoff law at node i and di denotes the neighborhood of i. This relation is exact for trees and can 
be considered as approximately correct for locally tree-like graphs, such as sparse random ones f37’, 139], where correlations 
between neighbors of a given node decay exponentially. One can develop further this recursion to obtain a set of three types of 
message-passing equations, one for each type of node, i.e. for each value of A' 1 . A fixed point of these equations can be found 
by iteration from arbitrary initial values for the messages until convergence. Then, one can collect at each edge the incoming 
and the outgoing converged messages to find the optimal configuration {7* ; .}(, y ) £ £ such that: 

I*j = arg minj {e u (T) + E y/ (-7) - /(||7||)} (6) 

where the last term is subtracted to avoid double counting of the cost of the single edge (ij). 


A. The mapping into a weighted matching problem. 

The min-sum algorithm as in 0 presents a computational bottleneck coming from the fact that for each output current 7 ly there 
is a large number of possible neighborhood’s configurations {fki\k&ii\j that are consistent both with the edge-disjoint constraints 
and with Kirchhoff’s law. In the calculation of the minimum in (0 one needs in fact to consider all possible combinations of 
paths entering and exiting node ;; the number of such combinations grows exponentially with the degree of node i. Nevertheless, 
the calculation can be performed efficiently by reducing it to a maximum weight matching problem 14611 on an auxiliary weighted 
complete graph G'. The nodes of G' are the neighbors k e di and the (symmetric) weights matrix Q will be defined as 

Qki = - min {E kl (v) + E u {-V)} + E ki ( 0) + £/,(()) (7) 
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FIG. 3: Mapping into a weighted matching problem. Left: intermediate step where G\ is built. On the leftmost part we show an 
example of several communications passing along (ij ) and exiting along the remaining neighbors k e di \ j. Right: the final 
step where G" 7 is built; the best configuration around node i when the blue current passes through (i j) is given by the minimum 
weighted matching on the complete auxiliary graph G" 7 . Edges red and green represent the best matching, i.e. the configuration 

where two other communications enter/exit neighbors of i \ j. 


where Eki(v) = Euihi) with s S MiV for v > 0, = —8^ v for v < 0 and = 0 for v = 0. Notice that this notation maps the 

M-dimensional vectors /, ; to the 2 M + 1 possible current configurations v allowed by the edge-disjointness constraint along a 
given edge. The computation of matrix Q, that requires 0(Mk 2 ) operations, should be performed only once at the beginning of 
the update routine for node i e G. 

Consider now a neighbor j e di and a given /r passing through edge (ij), we want to update £j 7 (/, 7 ). Assuming to know the 
other vertex 1 e di\j where the current /; entering (resp. exiting) node i can exit (resp. enter), then the least costly configuration 
in the remaining neighborhood is given by 

qj R = -M jl+ Yj (8) 

kedi\U,l) 

where Mji is the maximum weight of a matching on a complete graph G" 7 with k — 2 nodes, built from G' by removing nodes 
j and / (and all their incident edges). Recall that a matching is a subset of edges of G" 7 that do not share any vertex ||46j]. This is 
indeed equivalent to assigning to some of the remaining pairs of neighboring nodes currents v e [ -M ,..., M] that enters through 
one of them and exits through the other, such that the overall cost of the configuration is minimum. The key point is that the 
matching condition, i.e. the fact that edges in the solution set cannot have a vertex in common, in our problem translates in 
the condition of forbidding edge overlaps. Hence, thanks to this auxiliary mapping, we are able to reduce the computation of 
the update rule for the MP equations of the edge disjoint path problem to the solution of a standard (polynomial) combinatorial 
optimization problem, i.e. maximum weight matching. In Figure 0 we give a diagrammatic representation of the mapping. 
Note that in the maximum matching problem, edges in the input graph with negative weight can be simply removed. Notice 
that the neighboring current v can also be a priori equal to /r in this algorithm, because the configurations where jj appears in 
more than one pair of edges will be eliminated in the minimization calculation as they have higher cost in our formulation. The 
minimum weight is thus independent of /j, i.e. of which message we are updating, a fact that allows reducing the complexity of 
the algorithm by a factor M. 

Finally one needs to minimize over / given the matrix q mm : 

E'j 1 (^) = rnn [e' h (^ + qj"\ + c u (ju) (9) 

where c,- 7 (ju) is the cost of edge (i j), that in our case is 0 if/./ = 0 and 1 otherwise. We can notice that in order to evaluate each 
term inside the brackets we need to perform a matching optimization on each of the (k - 2)-node complete graphs G" 7 built 
V / e di \j. Each of these matching routine has complexity 0(k? log k) [47J and there are 0(k 2 ) possible combinations of j and 
l. Reminding that we first need to evaluate the weight matrix Q, the overall complexity of this algorithm will be: 

0(k 5 log k + Mk 2 ) 

which is polynomial in the variables k and M. Once we have performed this whole procedure, we get all the information we 
need to calculate the 2 M + 1 update messages E’k '(//), for each j e d i, adding a term O(kM) to the final complexity (which is 
nonetheless negligible compared to the previous two). 

The case of yU = 0, in which no current passes through edge (if) regardless of what happens on the other edges, is addressed 
by calculating a matching on the (k - 1 )-node complete graph composed of all nodes l 6 d i\j. If i is either a sender or a receiver, 
i.e. A. 1 ' 6 (±1) for a given fi e [1,..., M\, the same computation can be performed provided that an auxiliary node, indexed by 
the communication label /./ is added to the original graph G and connected to node i, such that its exiting messages will be fixed 
once at the beginning in the following way and never updated: E^fv) = —oo if 0 < v = n (sender) or 0 < —v - fi (receiver), and 
E,jj(v) = +oo otherwise. 
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B. The role of reinforcement. 

In order to aid and speed-up convergence of the MP equations, we used a reinforcement technique |48l|49}], in which a set 
of external local fields h\-(p) = Ek{p) + E’.^-p) - c,j(p) act on the messages gradually biasing them to align with themselves. 
The reinforcement is introduced by promoting edge costs to become communication-dependent quantities defined as linear 
combinations of the cost at the previous time-step and the reinforcement local fields: 

Cy V) = CyO“) + yth\jUi) (10) 

with c?- {fi) = Cij. This cost will then be inserted into equation © to replace the term Cjj(p). This has the effect to lead the 
messages to converge faster, gradually bootstrapping the system into a simpler one with large external fields. In practice we 
choose y, = ip and one has to choose the growth rate of y by tuning the reinforcement parameter p, that controls the trade-off 
between having a faster convergence and reaching a better solution. We tested p on instances on three types of graphs to finally 
choose to fix it to p = 0.002 in the rest of the simulations. In Figure[4]we could notice that this value achieves comparable results 
(inset) in terms of M acc /M to lower p in less time. 

In Figure [5jleft) we report the number of converged instances (over 100 realizations) for standard MP (without reinforcement) 
on four types of random graphs (as described in the next section) and fixed size V - 1000 and average degree ( k) - 3. The 
convergence failure of the standard MP increases considerably with M/V until it reaches a peak value, then it decreases. On the 
contrary, when reinforcement is used, convergence is always achieved in less than 100 steps (right panel). 
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FIG. 4: Reinforcement performance. Number of iterations to reach convergence as a function of reinforcement parameter p on 
BRITE graphs AS-BA217 (V=100) with M = 25,40, blrandl (V=500) with M = 125 and mesh 15x15 with M = 22. Inset: the 
number of accomodated paths M acc is substantially unchanged in the range of parameter values under study. 


IV. RESULTS ON RANDOM GRAPHS 

First, we tested the MP algorithm on various types of random graphs, with fixed size V = \'V\ = 1000 and average degree 
( k) = 3,5,7: regular random graphs (Reg), Erdos-Renyi random graphs (ER) |50'|. random graphs with power-law distribution 
(SF)il and a set of graphs (RER) obtained adding edges independently with probability p starting from a ko-regular random 
graph (for large V, the final average degree of such graphs is (k) = ko + d, with d = pV). We compared the performance 
with a multi-start greedy algorithm (MSG) |3^]. This heuristic algorithm calculates paths by iteratively choosing a (random) 
communication p, finding the corresponding shortest path and removing the edges belonging to the path from the graph. The 
process is repeated until either there are no paths left to be routed or no communications can be accommodated anymore in 
the graph. The multi-start version repeats the same procedure a given number of times and keeps the best solution in terms 
of M acc , the number of accommodated paths. A bounded-length version |52j] of MSG has been used to develop an iterative 
algorithm to solve the RWA using EDP in |fl8ll : its performance was comparable to the one obtained using a linear programming 
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FIG. 5: Left: Fraction of instances in which convergence standard MP fails (reinforced MP always converged in our 
experiments). Right: number of iterations for convergence for standard MP and reinforced MP (p = 0.002) in case of random 
graphs of V = 1000 and ( k) - 3 as a function of M/V. Notice how the reinforcement term, besides ensuring convergence, 

greatly improves the convergence time. 


solver on graphs of small sizes (V < 40) but with faster execution times. This makes it suited to be tested on larger graphs. A 
disadvantage of the greedy method is that it relies heavily on the order in which communications are accommodated (it disregards 
the information about sender-receiver pairs other than the ones already accommodated). The difference in the performances 
of the message-passing and greedy algorithm could then be used to assess the relevance of local information usage in such 
optimization problem. We tested both the standard multi-start and the bounded-length version but we found equal results with 
the first being slightly faster, in our tests, in terms of execution times. Thus we decided to use the standard MSG in our 
simulations. First we compared the results in terms of number of accommodated paths M acc by calculating the ratio M acc /M. In 
Figure[6]we show the behavior of M acc /M for each type of random graph and V = 10 3 , (k) = 3 using MP, reinforced MP and 
MSG. Both MP versions perform better than MSG, with the standard MP giving better results. The corresponding results for 
( k ) = 5 are similar (not reported) but the value M acc /M < 1 is reached at higher values of M/V and standard MP and MP with 
reinforcement give almost always the same solutions. The case ( k ) = 7 is not reported because, given the high number of edges, 
the solutions are often trivial (i.e. M acc /M = 1), a part from the case of SF graphs where we have instead M acc /M < 1 due to the 
presence of many small degree nodes. We also studied the total path length as a function of M/V for the solutions, obtained with 
the different algorithms. We consider the ratio between the total path lengths obtained with greedy and MP for solutions in which 
the number M acc of accommodated path is the same. In Figure [7] we can see that MP always outperforms the MSG algorithm 
for all types of graph under study. The results for the SF graph with (k) = 7 are quite different from the other graphs: both for 
MP and MSG the ratio departs from 1 at rather small values of M/V, possibly because the maximum number of accommodated 
paths is limited by the existence of many small degree nodes that act as bottlenecks, preventing the use of many alternative edge- 
disjoint routes. The scaling behavior of the fraction 1 - M acc /M of unaccommodated communications and the average total path 
length L/V of accommodated paths with the system size in the solutions obtained using the MP algorithm is shown in Figured 
for regular random graphs and ER random graphs. These quantities are plot as functions of the scaling variable x - — ’° E - . Note 
that when paths do not interact, x is a measure the total path length per site, as the average path lenght is proportional to log V. 
In the top panels, two regimes are visible: for small x, all communications can be accommodated, whereas at some value x* the 
curves for different values of V depart from zero. This behavior can be interpreted as a SAT/UNSAT transition, in analogy with 
the terminology of constraint-satisfaction problems EJ. The collapse of the curves L/V for different values of V is very good in 
the region in which all paths can be accommodated. On the contrary, in the UNSAT region, the curves for different sizes do not 
collapse anymore, though the relative difference between them seems to decrease by increasing the system size, and the curves 
for the largest graphs analyzed (V = 8000, 10000) are almost superimposed. We argue that x is the correct scaling variable in 
the limit of infinitely large graphs, and the observed mismatch could be due to finite-size effects. The change of slope in the 
roughly linear behavior of the average total length L/V is motivated by the fact that in the SAT region, all communications can 
be accommodated at the cost of taking longer paths with respect to those actually accommodated in the UNSAT region. 
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FIG. 6: MP vs greedy performance. We plot the performance in terms of M acc /M for (from top to bottom) regular, RER, ER 
and SF graphs of fixed size V = 10 3 and average degree (k) = 3. Error bars are smaller than the size of the symbols. 


V. COMPARISON WITH OTHER METHODS 


A comparison between the performances of the MP algorithm and those of alternative algorithms proposed in the literature 
Enm is reported in Table Q] As benchmark instances we used: two internet-like topologies generated using the BRITE graph 
generator [53] with parameters set as in (36} ]: mesh graphs of sizes 15x15 and 25x25, Steiner and planar graphs as reported 
in m. For each of these graphs we used the same set of sender-receiver pairs of size M = 0.10V, 0.25V, 0.40V used in l( 3lll . 
For each of these instances we ran the MP, MP with reinforcement and MSG algorithms 20 times and collected the average, 
minimum and maximum number of accommodated paths M acc along with the average computational time in seconds. All results 
are reported in Table Q] 


A. Other optimization methods. 

A part from the multi-start greedy, we used as comparison two more structured algorithms. The first one is an Ant Colony 
Optimization metaheuristic IBol . This method builds an EDP solution incrementally from partial solutions provided by a set 
of M ants. Each ant generates a path for a given communication making probabilistic decisions during the construction steps. 
These are made by processing local information modeled as pheromone information provided by other ants. The advantage 
of this method is to divide the EDP in subproblems and to use local information. The drawback is that it relies on several 
parameters that need to be carefully tuned in order to have a sensitive solution. Moreover the computational time increases 
considerably with the system size. The second algorithm is a Montecarlo-based Local Search m, that uses as main Montecarlo 
step a path rewiring based on rooted spanning trees. Unfortunately the running time grows rapidly with the system size, making 
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M/V 


M/V 


FIG. 7: Length performance. We plot (left) the relative performance of MSG over MP in terms of total length of the solution 
paths: y = 100 {L g /LMP - 1). Here L g and Lmp denote the total path lengths calculated with MSG and MP respectively. We use 
Reg, RER, ER and SF graphs of fixed size V = 10 3 and average degree (k) = 3,5,7 (from top to bottom). On the right we 
report the number of instances where the two algorithms find the same solution in term of M acc /M over 100 realizations. 


it computationally expensive when used on large graphs. Results are reported in Table Q] Finally, we performed simulation using 
the multi-start greedy heuristic described above. 


B. Results. 

In Table H] we report the performance comparison in terms of M acc between the two versions of MP (with and without rein¬ 
forcement) and the other 3 types of algorithms. The message-passing always performs equal or better than the other methods. 
Surprisingly the best performances are given for meshes and planar graphs, where we would expect the failure of MP due to 
the existence of short loops. What we find instead is that, even though the standard MP converges in few of these instances 
on meshes, the version with reinforcement always finds a solution that is always better than the other algorithms. The larger 
performance gap is seen on larger set of commodities and bigger graphs. Performance improvement reaches 27% with respect 
to LS, the best one between the other algorithms tested. The same considerations can be made in the case of planar graphs. We 
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FIG. 8: Finite-size effects. We plot 1 - M acc /M (top) and the total length per node L/V (bottom) for Reg (left) snd ER (right) 
graphs as a function of the scaling variable M1 ° gV . We can notice the finite-size effects decreasing with system size leading to 
the curves corresponding to the biggest graphs V = 8000, 10000 to almost superimpose. Note that in the SAT phase the total 
length grows linearly in log V for all system sizes as expected but in the UNSAT phase the graphs split. Error bars are smaller 

than point size. 


claim that this gap would increase with system size, but unfortunately the size of benchmark graphs remains limited to V < 500. 
Moreover these alternative algorithms do not consider path length optimization, thus we cannot compare the performance with 
respect to this variable. The ACO has been recently tested on several types of graphs (still with V < 500) against a Genetic 
Algorithm (GA) in [54]. It performed better than GA in the case of BRITE graphs 1-6 and 14% worse in the case of 10x10 
and 15x15 mesh graphs. The MP algorithm always outperforms ACO and in the case of 15x15 mesh the gap reaches 23.5%. 
Unfortunately neither the GA has been tested on larger graphs nor gives results in terms of path length of the solutions. 


VI. CONCLUSIONS 

The EDP problem is a combinatorial optimization problem that finds applications in several traffic engineering problems, 
from VLSI design to routing and access control management in communication networks. In this work we proposed a min- 
sum message-passing algorithm to find the maximum number of communications M acc that can be accommodated in a network 
subject to edge-disjoint constraints and minimizing total path length at the same time. We devised an efficient method to 
implement these equations by exploiting a mapping into a minimum weight matching problem on an auxiliary graph. The 
standard MP algorithm and the version with reinforcement consistently outperform alternative algorithms found in the literature 
on different types of benchmark graphs in terms of the fraction M acc /M of accommodated communications. We found two 
different behaviors: on some “easy” instances, all algorithm accommodate all requests, providing the same results and suggesting 
that these could be the optimal ones; there are non-trivial instances in which M acc /M < 1, but the message-passing algorithm 
always outperforms the other algorithms in terms of the number of accommodated paths. In particular we obtained better results 
in the case of meshes and planar graphs, even though these topologies are not locally tree-like as required by the cavity method. 
In these cases, we could always ensure convergence of the MP equations by exploiting a reinforcement technique. The quality 
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of solutions improves with decreasing the reinforcement parameter, such that we could always find better solutions than those 
obtained using the other algorithms under study. Unfortunately, for the heuristic algorithms employed on the benchmarks we 
could not access other relevant metrics such as the average total path length, as it was not considered before in the literature 
JUS. Nonetheless we could directly compute such quantity for a multi-start greedy heuristic in several graphs, finding that 
MP always gives a lower average path length for solutions with the same fraction of accommodated communications. 

In conclusion, combining the good performance results, in terms of traffic and path length, with the polynomial time imple¬ 
mentation, the use of the MP algorithm opens new perspectives in the solution of relevant routing problems over communication 
networks such as the RWA in optical networks. In particular, it would be interesting to apply the MP algorithm in the iterative 
construction of RWA solutions over communication networks with finite link capacity, as it has been done for other types of 
EDP algorithms. 
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Appendix A: Convergence criterion 

Given a decision variable d' to be calculated at each iteration update step t, an integer variable n and a time step T max we have 
convergence if, for n consecutive iteration steps, d' does not change, and we fix a maximum iteration time T max to update MP 
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equations. Formally this writes: 

3t 0 e[l,T max -n] s.t. d' 0+i = d'° Vi=l,...,n (Al) 

In our simulations we defined the decision variable as the total difference of the optimal currents (calculated edge by edge) 
between two consecutive iteration steps : 


d ' = 2 

ftfleE 


- 


(A2) 


where /rF = | min /J= _ M M {^(/r) + £_/,■(—//) - c i; (/r)J | and convergence is reached when d‘ - 0 for n consecutive time steps. 


Appendix B: Benchmark results 
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285 

5.70 

37.00 

37 

37 

37.05 

37 

38 

35.80 

35 

37 

34.00 

33 

36 

36.00 

35 

37 

2.70 

5.56 

2.70 

plan200 10 

200 

583 

5.83 

20.00 

20 

20 

20.00 

20 

20 

20.00 

20 

20 

20.00 

20 

20 

20.00 

20 

20 

0.00 

0.00 

0.00 

pan200 25 

200 

583 

5.83 

- 

48.95 

48 

50 

46.50 

46 

48 

41.80 

39 

43 

45.95 

45 

48 

4.17 

16.28 

4.17 

plan200 40 

200 

583 

5.83 

. 

60.65 

58 

62 

52.95 

52 

56 

49.35 

47 

51 

55.70 

54 

58 

10.71 

21.57 

6.90 

plan500 10 

500 

1477 

5.91 

50.00 

50 

50 

50.00 

50 

50 

50.00 

50 

50 

44.95 

42 

47 

50.00 

50 

50 

0.00 

6.38 

0.00 

plan500 25 

500 

1477 

5.91 

. 

92.29 

90 

94 

78.15 

76 

80 

60.95 

57 

65 

78.20 

77 

80 

17.50 

44.62 

17.50 

plan500 40 

500 

1477 

5.91 

- 

122.31 

119 

124 

92.60 

90 

95 

82.85 

78 

86 

100.15 

97 

102 

30.53 

44.19 

21.57 


TABLE I: Message-passing and multi-start greedy performances. Columns 1-4 give the characteristics of the benchmark. For 
each algorithm, columns 1-3 represent the average, the minimum and the max number of accommodated paths over 20 runs of 
a given set of commodity instance respectively. ACO and LS performances are reported in lf3ll[3d l. Performance comparison 
between MP and the other algorithms is given in the three last columns, representing the performance ratio 
100 ■ (M^/Macc - 1) where alg indicates the algorithm used (MSG, ACO and LS respectively). We use as the best one 

between MP with and without reinforcement. 































































